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Abstract 

Binary Bargmann symmetry constraints are applied to decompose soliton equa- 
tions into finite-dimensional Liouville integrable Hamiltonian systems, generated 
from so-called constrained flows. The resulting constraints on the potentials of soli- 
ton equations give rise to involutive solutions to soliton equations, and thus the in- 
tegrability by quadratures are shown for soliton equations by the constrained flows. 
The multi-wave interaction equations associated with the 3x3 matrix AKNS spec- 
tral problem are chosen as an illustrative example to carry out binary Bargmann 
symmetry constraints. The Lax representations and the corresponding r-matrix 
formulation are established for the constrained flows of the multi-wave interaction 
equations, and the integrals of motion generated from the Lax representations are 
utilized to show the Liouville integrability for the resulting constrained flows. Fi- 
nally, involutive solutions to the multi-wave interaction equations are presented. 
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1 Introduction 



The nonlinearization technique has been attracting broad attention in the 
soliton community [1-3] . It provides a powerful approach for analyzing soliton 
equations, in both continuous and discrete cases, especially for showing the 
integrability by quadratures for soliton equations. A kind of specific Bargmann 
symmetry constraints [4] yields the required constraints on the potentials of 
soliton equations in nonlinearization. One of two conserved covariants involved 
in each Bargmann symmetry constraint is Lie point, but the other is not Lie 
point, contact or Lie-Backlund, which is expressed in terms of eigenf unctions 
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and adjoint eigenfunctions of spectral problems associated with soliton equa- 
tions. Such symmetry constraints nonlinearize spectral problems of soliton 
equations into finite-dimensional Liouville integrable systems, and the con- 
straints on the potentials present involutive solutions of soliton equations. 

Binary Bargmann symmetry constraints also establish a bridge between 
infinite-dimensional soliton equations and finite-dimensional Liouville inte- 
grable systems [4]. The constrained flows resulted from spectral problems of 
soliton equations are utilized to show the integrability by quadratures for soli- 
ton equations [5,6] and pave a method of separation of variables for soliton 
equations [7,8]. The study of binary symmetry constraints itself is an impor- 
tant part of the kernel of the mathematical theory of nonlinearization. 

This paper is structured as follows. In section 2, a general theory is recalled 
for reference. Then starting from section 3, an illustrative example is carried 
out to manipulate binary Bargmann symmetry constraints. Finally in section 
8, a summary is given along with concluding remarks. 



2 General theory 

Let us address a skeleton theory of binary Bargmann symmetry constraints 
(for example, see [9] for details). We start from a matrix spectral problem 

<j> x = U<f> = U{U, A)0, U = (E^-)rxr, = {(f),, • • • , r ) T (2.1) 

with a spectral parameter A and a potential u — (ui, • • • , u q ) T . Suppose that 
the compatability conditions U tm — + [U, V^] = 0, m > 0, of the spectral 
problem (2.1) and the associated spectral problems 

tm = ^ (m V = ^ {m) K^,---;A)0, vw = (yfpUr, m>o, (2.2) 

determine an isospectral (X tm = 0) soliton hierarchy 

u tm = X m {u) = JG m = J-^, m > 0, (2.3) 

where the Hamiltonian operator J and the Hamiltonian functionals H m can 
be determined by a trace identity. Evidently, the compatability conditions of 
the adjoint spectral problem and the adjoint associated spectral problems 

^ = -f/ T K Vt m = -V^ T (u, X)tp, m > 0, (2.4) 
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where ip — (ipi, • ■ • , -0v) T , still give rise to the same soliton hierarchy (2.3). 

t i i i r.i i r-i-1 t ,t9U(u,X) , . 

It has been pointed out 4 that J— = E Jip is a common 

ou ou 
symmetry of all equations in the hierarchy (2.3), where E is the normalized 

constant. Upon introducing iV distinct eigenvalues Ai, • • • , Ajv, we obtain 

0^ = U(u, X s )^ s \ ^ = -U T (u, A S )V (S) , 1 < s < N; (2.5) 
4 S 2 = V^ m \u, A s )0«, 4 S 2 = -V im)T (u, X s )^ s \ l<s<N, (2.6) 

where the corresponding eigenfunctions and adjoint eigenfunctions are denoted 
by 0^ and ip^ s \ 1 < s < N.It is assumed that the conserved covariant G mo is 
Lie point, and then the so-called binary Bargmann symmetry constraint reads 

X rno = EEsJ 5 ^, i.e., JG mo = JX> (S)T ^|^ (S) , (2.7) 

where E s , 1 < s < N, are N normalized constants. The right-hand side 
of the symmetry constraint (2.7) is a linear combination of iV symmetries 
E s JSX s /Su, 1 < s < N. Such symmetries (or the corresponding conserved 
covariants 5X s /5u, 1 < s < N) are not Lie point, contact or Lie-Backlund 
symmetries, since 0^ and ip^ can not be expressed in terms of x, u and 
derivatives of u with respect to x to some finite order. 

Suppose that the symmetry constraint (2.7) has an inverse function 

M = S(0«,...,0W;^),...,^)). (2.8) 

Replacing u with u in the system (2.5) or the system (2.6), we obtain the 
so-called spatial constrained flow and the so-called temporal constrained flow: 

<f>^ = U(u, X s )^ s \ = ~U T (u, A S )V (S) , 1 < s < N; (2.9) 

<$> = V™(u, A s )0« 4i = -V^ T (u, X s )^ s \ l<s<N. (2.10) 

The former is a system of ordinary differential equations, and the latter is 
usually a system of partial differential equations but it can be transformed 
into a system of ordinary differential equations under the control of (2.9). 

The main question of symmetry constraints is to show that the spatial con- 
strained flow (2.9) and the temporal constrained flows (2.10) under the control 
of (2.9) are Liouville integrable. Then if <f>^ and ip( s \ 1 < s < N, solve two 
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constrained flows (2.9) and (2.10) simultaneously, u = u will present a solu- 
tion to the mth soliton equation u tm = X m (u). It also follows that the soliton 
equation u tjn = X m {u) is decomposed into two finite-dimensional Liouville 
integrable systems, and the constraint u = u gives a Backlund transforma- 
tion between soliton equations and the resulting finite-dimensional Liouville 
integrable systems [10]. The whole process to carry out binary Bargmann 
symmetry constraints is called binary nonlinearization [11,5]. 



3 Multiwave interaction equations 



Let us begin with the 3x3 matrix AKNS spectral problem: 



<l> x = U{u,\)<l>, U 



( \ \ 

Ot\\ U 12 Ui 3 

u 2 i a 2 \ u 23 
U 3 i u 32 "3^ J 



<P2 



, (3.1) 



where U = diag(ai, a 2 , a 3 ), and oti,a 2 ,a 3 are distinct constants, and the 
potential u is defined by 



u = p(Ui) = (u 21 ,u 12 ,u 31 ,u 13 ,u 32 ,u 23 ) J 



(3.2) 



We introduce an associated spectral problem: 



PiA v 12 u 13 



<i> t = V(u,\)<j>, V 



v 21 (3 2 \ v 23 
v 3 i v 32 p 3 \ 



W + V u 



where V = diag(/3i, f3 2 , /3 3 ), and Pi, fa, fa are distinct constants. Under the 
isospectral condition X t = 0, the compatibility condition U t — V x + [U, V] = 
is equivalent to 



U u - Vi x + [U u Vi\ = 0, [U , Vi\ = [V , Ui\. 



(3.3) 



They give rise to the multi-wave interaction equations 



Uij, t = 



fa-P 3 



-u 



ij.x 



+ E 



3 Pk-fa Pk-Pj 



)u ik u kj , 1 < i ^ j < 3, (3.4) 
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which contain three-wave interaction equations arising in fluid dynamics and 
plasma physics [12,13], with U\ being chosen to be an anti-Hermitian matrix. 
Note that the compatability condition of the adjoint spectral problem and 
adjoint associated spectral problem: 

iP x = -U T (u,\)ij, A = -V T (u,\)ij, iP = (Vi,^3) T (3.5) 



still gives rise to the above multi-wave interaction equations. It is easy to find 
that the multi-wave interaction equations (3.4) have the Hamiltonian structure 

5H 

Uij,t = J 1 < i ^ j < 3, (3.6) 



where the Hamiltonian operator J and the Hamiltonian H are defined by 



J = diag((ai - a 2 )a , («i - a 3 )a , (a 2 - a 3 )cro), a 



■1 



, (3.7) 



H = f H dx, H = — — -(— — — - — — —)(u 12 u 23 u 31 



+U 2 iU 13 U 32 ) + 



3(«i — a 2 ) v «3 — a 2 a 3 — a\ 
Pi -fa 



2(«i - a 2 ) : 



(ui 2 u 21:X - u 21 u 12:X ) + cycle(l, 2, 3). (3.* 



4 Binary Bargmann symmetry constraints 

Let us consider binary Lax systems of the multi-wave interaction equations 
(3.4): 

<P X = U(u, A)0, i> x = -U T (u, (4.1) 
(j>t = V(u, A)0, i/h = -V T (u, (4.2) 



Based on the Hamiltonian theory, the conserved functional X(u) gives rise to 
a symmetry of (3.4): 

EJ 6 ^ = Ep([U ,p~\^)])=p([U ,p^V^^m, (4-3) 
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where E is the normalized constant and p is the mapping defined by (3.2). 
Having introduced N distinct eigenvalues Ai, • • • , Ajv, we obtain 

= U(u, A.)tf« ^ = -U T (u, X s )^ s \ l<s<N; (4.4) 
4 S) = V(u, \ 9 )<j>M, 4 S) = ~V T (u, A s )^ (s) , 1 < s < N, (4.5) 

where 0^ and ip^ are assumed to be denoted by 

= (01s, 02s, 03sf, i> (8) = ^2s, ^3sf, 1<S<N. (4.6) 

This multi-eigenvalue case yields a general symmetry of (3.4) (not Lie point, 
contact or Lie Backlund type): 

Zo :=p([a„,^(E^^#)])=p(Po,EA (s)r ]). (4.7) 

s=l OU s=l 

On the other hand, it is easy to find that the multi-wave interaction equa- 
tions (3.4) have a Lie point symmetry 

Y := p([T, U x \), V = diag( 7 i,72,73), 7* £lj,l<i±3< 3, (4.8) 

where 71,72,73 are constants. Actually, we can directly prove that Y and Z 
are two symmetries of the multi-wave interaction equations (3.4). Set 

5U X = [r, U x \ or [U , J2 (S V (S)T ], (4.9) 

s=l 

then [t/o,5Vi] = [Vb,5£/i] determines 

N 

5V 1 = [r, Vi] or [V , (S V (S)T ]- (4-10) 

s=l 

Now in view of (3.3), the symmetry problem requires to show that 

N N 

m,Wi) = ([^c/iljr,^]) or ([c/o,E^V (s)T ],[H,E^ (s V (s)T ]) 

s=l s=l 

satisfies the linearized system of the multi-wave interaction equations (3.4): 
(SUjt - (5V,), + [SU^Vi] + [U 1 ,5V 1 ] = 0, (4.11) 
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which just needs a direct computation. 

Therefore, binary Bargmann symmetry constraint of (3.4) reads as 



N 



Y = Z , i.e., [r,fA] = [[/ ,E^ (S V (S)T ]- (4.12) 



s=l 



Solving this for u, we have the constraints on the potentials 



u = m ;= — l(<f> h # l < j < 3 , (4.13) 

7i - 7j 



where (•, •) denotes the standard inner-product of ~R N , and 

= (0 a , • • • , <fev) T , tf< = (^i, • • • , ?M T , 1 < i < 3. (4.14) 

The substitution of w with u = (u^i, Uu, U31, U13, U23, ^32) in binary Lax 
systems (4.4) and (4.5) leads to the so-called constrained flows: 

= U(u, X s )^ s \ = -f/ T («, A s )^ (s) , 1 < a < N; (4.15) 
# = A a )0« ^ t (s) = -V^ t (m, A,)^ (s) , 1 < s < N. (4.16) 



These are two systems of ordinary differential equations. In the following two 
sections, we would like to show their Liouville integrability. 



5 Hamiltonian structures and Lax representations 



In order to prove that two constrained flows are Liouville integrable, we 
need generate their Hamiltonian structures and integrals of motion. First, the 
spatial constrained flow (4.15) can be easily expressed as the Hamiltonian form 

dH x dH x 



with the Hamiltonian 

H* = -j2«k{A* k ,* k )- £ ^r^<*fc, *i><*i, (5-2) 

fc=l l<fc</<3 " k P l 
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where A = diag(Ai, A2, • • • , Xn)- Second, the temporal constrained flow (4.16) 
has a similar Hamiltonian structure 



&it. = — — , = 1 < i < 3, 



(5.3) 



with the Hamiltonian 



fe=l l<fe</<3 7fc 7Z 



(5.4) 



It is known that Lax representations are basic objects to generate integrals of 
motion. In order to present Lax representations of the constrained flows (4.15) 
and (4.16), let us define a Lax operator L(X) by 



L(X) 



7i + E 



N M 



N foii 



E 



N 



1 X — Xi l=1 X — Xi 

N 



1=1 x - x l 



\ 



^ A- A, 

<Mll 

V SA-A, 



N 

E 



4>2li>3l 



1 X — Xi i= i x — Xi 



N 

E 



21 



IX- Xt 



N 



73 + E 



fol^3l 



\ X - Xi J 



(5.5) 



Now if we set U(X) = U(u,X) and V"(A) = V(u,X), then it is easy to check 
the following Lax representations: 

(L(A)) X = [U(X), L(A)], (L(X)) t = [V(X), L(X)}, (5.6) 



for the spatial constrained flow (4.15) and the temporal constrained flow 
(4.16), respectively. These Lax representations will be applied to construct the 
required integrals of motion for the Liouville integrability of two constrained 
flows in the next section. 



6 r-matrix formulation and Liouville integrability 



Let us specify the Poisson bracket: 

f '* e ^>- (6 - l) 
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An r-matrix formulation can be directly shown for the Lax operator L(X) 
defined by (5.5). 

Theorem 6.1 The Lax operator L(X) given by (5.5) has the following r- 
matrix formulation 



{L(X) f L(fx)} — 



1 



fi — X 



V,L 1 (X) + L 2 ( f i) 



3 



,V=Y,Eij® Eji, (6.2) 



where L±(X) = L(X) ® J 3 , L 2 (X) — I 3 <g) L(fi), and 

(Eij)ki = S ik 6ji, {L(X) ® L(n)} ijM = {(L(A)) ifc , 1 < i,j,k,l < 3. 

First from the Lax representations in (5.6), we have 
(ul 3 - L(X)) X = [U(X), vh - L(X)}, {vh - L{X)) t = [V(X), ul 3 - L(X)}, 

where v is a constant, and thus [14] 

(det(z// 3 - L(X))) X = 0, (det(z// 3 - L(X))) t = 0. 

Second from the r-matrix formulation (6.2), we have [15] 

{trL k (X),trL l (/2)} = 0, k,l>l. 

Expand the determinant of the matrix vl 3 — L(X) as 

det(z// 3 - L(X)) = v 3 - ^(Xy + F 2 {X)v - JF 3 (A), (6.3) 

where by Newton's identities on elementary symmetric polynomials, 
^i(A) = trL(X), F 2 (X) = \{{trL{X) f - trL 2 (X)), 
F 3 (X) = l(trL(X)) 3 + |trL 3 (A) - \{trL{X))trL\X). 

Thus, it follows that 

(J r i(A)) x = 0, (^(A)) t = 0, and {^(A),^)} = 0, 1 < i, j < 3. (6.4) 

Further expand J~i{X) as 

Ti{X) =Y,F a \- 1 , 1<2<3, (6.5) 

l>0 
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where obviously the F i0 's are constants. Then from the last equality in (6.4), 
we have 

{Fu, F jk } = 0, 1 < i,j < 3, k,l > 0. (6.6) 

The first two equalities in (6.4) and the equality (6.6) implies that two con- 
strained flows (4.15) and (4.16) have the common involutive integrals of mo- 
tion: Fu, 1 < % < 3, I > 1. Now it is a direct computation to verify the 
following theorem on the Liouville integrability of the constrained flows (4.15) 
and (4.16). 

Theorem 6.2 Two constrained flows (4-15) and (4-16) have the common in- 
volutive integrals of motion: Fu, 1 < % < 3, I > 1, of which the functions 
Fi S , 1 < i < 3, 1 < s < N, are functionally independent over a dense open 
subset of~R 6N . Therefore, the constrained flows (4-15) and (4-16) are Liouville 
integrable. 



7 Involutive solutions 

Since under the constraints on the potentials (4.13), the compatability con- 
dition of (4.4) and (4.5) is still the multi-wave interaction equations (3.4), 
solutions ($j(x, t), \l/j(x, t)) to the constrained flows (4.15) and (4.16) present 
solutions to the multi-wave interaction equations (3.4): 

Uij{x, t) = J -(^(x, t), t)), l<i^j<3. (7.1) 

1% ~ lj 

This also shows the integrability by quadratures for the multi-wave interac- 
tion equations (3.4) since and can be determined by quadratures. On 
the other hand, it is easily found that under the Poisson bracket (6.1), two 
Hamiltonians H x and if* commute, i.e., 

Hence, the above solutions given by (7.1) determine involutive solutions to the 
multi-wave interaction equations (3.4). If we denote two Hamiltonian flows of 
(4.15) and (4.16) by g^ x and gf respectively, then we have 

it ij 

= ^^(gfgT^gfgT^^), i < i± j < 3, 
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where the initial values <E>io and \J/jo of and can be taken as any arbitrary- 
constant vectors of M. N . 

Summing up, solutions of the constrained flows (4.15) and (4.16) lead to 
involutive solutions of the multi-wave interaction equations (3.4). Moreover, 
such involutive solutions show us (a) the richness of solutions of the multi- 
wave interaction equations (3.4), and (b) the integrability by quadratures for 
the multi-wave interaction equations (3.4). 



8 Summary and conclusion 



Binary Bargmann symmetry constraints in the continuous case decompose 
soliton equations (PDEs) into constrained flows (ODEs). The r- matrix formu- 
lation can be used to show the Liouville integrability of the constrained flows. 
The resulting constraints on the potentials give rise to involutive solutions 
to soliton equations and thus show the integrability by quadratures for soli- 
ton equations. The whole process of binary Bargmann symmetry constraints, 
called binary nonlinearization, can be depicted as follows. 



Lax system and 
adjoint Lax system 



<p x = U(u, \)<f> 

<t> t = v{u,\)4> 
i>x = -u T (u, 

^ t = -V T {u,\)i> 



Integrability by 
quadratures for 
U t -V x +[U,V]=0 



by 



symmetry 
constraints 



by 



involutive 
solutions 



Constrained spatial 
and temporal flows 



(j) x = U(u(<l>,ip),X)<f> 
0, = \/(w(0,^),A)0 



by r-matrix 
formulation 



finite-dimensional 
Liouville integrable 
Hamiltonian systems 



The multi-wave interaction equations (3.4) have been taken as an illus- 
trative example. Binary Bargmann symmetry constraints (4.12), containing a 
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set of arbitrary distinct constants 71, 72, 73, were proposed for the multi-wave 
interaction equations (3.4). Two finite-dimensional Liouville integrable Hamil- 
tonian systems (5.1) and (5.3) resulted from the constrained flows determine 
involutive solutions to the multi-wave interaction equations (3.4). Our result 
with a special case T = Vq, i.e., diag(7i, 72, 73) = diag(/3i, /3 2 , fa), gives rise to 
all the results established in [16]. 

Of special interest in the study of binary symmetry constraints are to cre- 
ate new classical integrable Hamiltonian systems which supplement the known 
classes of classical integrable systems [17] and to expose the integrability by 
quadratures for soliton equations by using constrained flows [9]. We point out 
that high-order symmetry constraints with involved symmetries having de- 
generate Hamiltonians need particular consideration [18]. A profound math- 
ematical theory on binary symmetry constraints, especially on Hamiltonian 
structures of constrained flows, will be discussed elsewhere. 
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